Two-dimensional Mott-Hubbard electrons in an artificial honeycomb lattice 
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Electrons in artificial lattices enable explorations of the impact of repulsive Coulomb interactions 
in a tunable system. We have trapped two-dimensional electrons belonging to a gallium arsenide 
quantum well in a nanofabricated lattice with honeycomb geometry. We probe the excitation spec- 
trum in a magnetic field identifying novel collective modes that emerge from the Coulomb interaction 
in the artificial lattice as predicted by the Mott-Hubbard model. These observations allow us to 
determine the Hubbard gap and suggest the existence of a novel Coulomb-driven ground state. This 
approach offers new venues for the study of quantum phenomena in a controllable solid-state system. 



When electrons roam in a solid they experience the 
crystalline potential created by a periodic arrange- 
ment of coupled quantum units (such as ions, atoms, 
or molecules). Numerous properties of solids can 
be explained within the paradigmatic theory of Bloch 
bands [1], which neglects Coulomb interactions. These, 
however, often lead to profound qualitative changes that 
are particularly pronounced in solids with narrow en- 
ergy bands [2-4]. These strongly correlated materials 
display exotic ordering phenomena and metal-insulator 
phase transitions [3, 4]. 

Mott showed [5] that interaction-induced insulators are 
better described in real space (rather than in momentum 
space), in which the solid is viewed as a collection of 
localized electrons bound to atoms with partially filled 
shells. Electrons hopping through the lattice are ab- 
sorbed and emitted from the atoms thus originating two 
bands, which are split by the energy cost of having two 
electrons with antiparallel spin on the same atomic site. 

Hubbard subsequently introduced a model Hamilto- 
nian with on-site interactions that displays split bands 
{Hubbard bands or HBs) in the strongly correlated (or 
atomic) limit [6]. HBs (and their coexistence with 
quasiparticle bands in the correlated metallic phase) are 
"quintessential" features of such strongly correlated sys- 
tems. Fermions on a honeycomb lattice, in particular, 
have been predicted to display unusual correlated phases 
of matter such as topological Mott insulating [7] and 
quantum spin liquid phases [8]. Thus the creation of 
artificial systems with a high degree of tunability that 
offer access to the Mott-Hubbard physics is an extremely 
appealing endeavor [9]. 

Here we report on the creation of an artificial lattice 
with honeycomb geometry for electrons, and we demon- 
strate the formation of HBs due to strong correlations. 
The artificial lattice is obtained by nanofabrication ap- 
plied to a GaAs heterostructure that supports a high 
quality two-dimensional electron gas (2DEG) [10-12]. 



We probe the excitation spectrum of the electrons by in- 
elastic light scattering and observe signatures stemming 
from strong Coulomb interactions, which are tuned by 
the application of an external magnetic field. 

We find that carriers in the patterned structures sup- 
port a novel collective mode that has energy that scales 
like V^, where B is the component of the magnetic field 
perpendicular to the 2DEG. A theoretical analysis based 
on a minimal Hubbard model reveals that the mode en- 
ergy is determined by the on-site Coulomb interaction 
and thus represents direct evidence of the existence of 
HBs in the 2DEG subjected to the artificial lattice. 

At low temperatures we find evidence for the opening 
of an unexpected gap in the spin excitation spectrum at 
large B fields. This observation underpins the impact 
of collective phenomena in artificial lattices. We argue 
that the observed gap reveals the occurrence of a new 
correlated phase of electrons in a honeycomb lattice akin 
to one of those discussed in the context of graphene at 
high magnetic fields [13-15]. These findings pave the 
way for the possibility to explore graphene-like physics in 
the ultra-high-magnetic-field limit in which the magnetic 
length is smaller than the lattice constant of the artificial 
crystal - a regime not accessible in graphene. 

The capability of observing Mott-Hubbard physics in 
nanostructured semiconductor devices with honeycomb 
geometry may open new approaches for the investiga- 
tion of quantum phases of strongly correlated condensed- 
matter systems. Given that the interaction strengths 
governing the physics of the 2DEG can be finely tuned 
by design and the application of external electric and 
magnetic fields [16, 17], such scalable solid-state systems 
offer great promise to further expand the current realms 
of study offered by quantum emulators that have been so 
far realized with cold atom gases in optical lattices [18- 
20]. 

The sample used in this study is the host of a 
2DEG in a 25 nm wide, one-side modulation-doped 
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FIG. 1: a) Scanning electron microscopy (SEM) image of the semiconductor artificial lattice. An expanded view of the SEM 
image showing a single honeycomb cell (2r ^ 60 nm, a ^ 130 nm). The two-dimensional electron gas is positioned 170 nm 
below the surface with a low-temperature mobility of 2.7 x 10^ cm^/(Vs). We also sketch a cartoon of the two-dimensional 
potential trap for electrons induced by the nanofabricated pillar at the surface, b) Geometry of the light scattering experiment: 
C(;l,s labels the incident (scattered) photon energy and B = 5° is the tilt angle, c) Resonant inelastic light scattering spectra 
showing the cyclotron mode and the new low- lying collective mode at B = 5.48 T and T = 1.7 K. d) Evolution of the energies 
of the cyclotron mode (black filled circles) and of the new collective mode at frequencies cjhb (red filled squares) at T = 1.7 K. 
The black dashed line is a linear fit to the data using huuc = heB/{m*c). We find m* = 0.067 rrie with rrie the bare electron 
mass, in agreement with the bulk GaAs value. The red dashed line is a fit with Hujub = ay^BfTJ and a ^ 2 meV. 



Alo.iGao.gAs/GaAs quantum well. The procedures for 
nanofabricating the artificial lattice [10, 21] are de- 
tailed in the Supplementary Material (SM). The artifi- 
cial honeycomb lattice extends over a 100 /im x 100 /im 
square region with a lattice constant a ~ 130 nm [see 
Fig. la)]. We denote by Vq the amplitude of the arti- 
ficial lattice potential. Here we focus on a sample with 
an estimated [11] Vq ^ 4 meV and electron density after 
processing ne 3 — 4 x 10^^ cm~^ corresponding to an 
average number of eight electrons per site. 

The inelastic light scattering experiments were per- 
formed in a backscattering configuration [see Fig. lb)] in 
the temperature range 50 mK — 4 K. The light scatter- 
ing technique gives direct access to the collective modes 
of the system that manifest as sharp peaks in the inten- 
sity of the scattered light at a given energy shift from the 
laser energy. Crucial for these observations is the reso- 
nant enhancement of the light scattering cross section 
that occurs as the incident laser energy is scanned across 
an inter-band transition of the host GaAs semiconductor 
(see SM). 

The nanostructured 2DEGs displays well-resolved 
quantum Hall signatures below 3 T with the honey- 
comb potential manifesting itself in a modulation of the 
magneto-resistivity periodic in B (data not shown and 
Ref. 11). At larger fields, an increase in the longitudinal 
resistivity signals a crossover to a regime of suppressed 
inter-site hopping in which new collective modes emerge. 



In addition to the ordinary cyclotron mode [black filled 
circles in Fig. Id)] at energy = heB /{rrt'c)^ being 
the GaAs electron effective mass [22], the light scattering 
spectra display a new mode at lower energies as shown 
in Figs. Ic) and Id). This new excitation reveals its col- 
lective character in the sharpness and intensity of the 
light scattering peak [23]. The surprising sublinear de- 
pendence of the energy of the new mode on B is shown 
in Fig. Id) (red filled squares). 

We identify the sublinear collective mode with a Hub- 
hard mode^ i.e. an excitation across split HBs. In the sim- 
plest scenario, this excitation emerges within the single- 
band Hubbard model [3, 4, 6] that assumes a maximum 
concentration of two electrons per site. We proceed by 
first evaluating the Mott-Hubbard excitation gap as a 
function of B and then we demonstrate that it weakly 
depends on electron concentration consistently with ex- 
perimental data shown in Fig. 3b). Similar conclusions 
can be reached by employing multiband generalizations 
of the Hubbard model (see SM). 

The single-band Hubbard Hamiltonian encodes a com- 
petition between two energy scales: the kinetic energy 
t, which measures the overlap between electronic wave- 
functions on neighboring lattice sites, and the interac- 
tion energy [/, which measures the strength of the on-site 
Coulomb repulsion between two electrons: 

?i = ^ c\cj + So ^ + ^ ^ rii^rii^ • (1) 
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FIG. 2: a) A cartoon of the spectral function A{io) of the 
patterned/unpatterned 2DEG (red/black). The Landau level 
peaks at cj = cjc + 1/2) are split by on-site Coulomb in- 
teractions into Hubbard lower and upper peaks, which are 
separated by U ^ /Ib, where Ib = y^Kc/JeB) is the mag- 
netic length, b) The relevant electronic process which con- 
tributes to the Raman scattering cross section. The initial 
state is labeled by |1), the final state by |2), while the inter- 
mediate state with one hole and an extra electron is labeled 
by \n). The final excited state is separated from the ground 
state by the Hubbard charge gap [/, i.e. by the energy cost 
of having two antiparallel spin electrons on the same site. In 
the intermediate state we have also depicted the absorbed (at 
frequency cjl) and emitted (at frequency cjs) photons. The 
square wells denote two neighboring minima of the artificial- 
lattice potential. The core levels are not shown. The green 
areas denote valence-band electrons, which are assumed to be 
unaffected by the periodic modulation. 



Here q (q) creates (destroys) an electron at site i (the 
sum in the first term is over all pairs of nearest-neighbor 
sites), and hi = c\ci is the local number operator; 
denotes the energy of the single state which is available 
at each site i: this can be either empty, singly, or dou- 
bly occupied. In writing Eq. (1) we have neglected first- 
neighbor {i.e. inter-site) interactions. In the atomic, 
strongly correlated limit U ^ two split HBs emerge 
out of a single narrow band [6]. More precisely, this 
means that for U ^ t the spectral function, i.e. the 



tunneling density-of- states, A{uj) of the model described 
by Eq. (1) develops two peaks, one at huo = Sq and one 
Sit huj = So -\- U . The emergence of HBs when the ratio 
U/t increases from the weakly to the strongly correlated 
regime is accurately described by dynamical mean-field 
theory [2-4]. 

In the experiments, the strongly correlated regime 
U t is achieved when B quenches the hopping am- 
plitude t and increases the interaction energy U. The 
Hubbard- /7 interaction scale can be written in terms of 
localized Wannier functions (j){r) as 

U = j d?r j dV'|0(r)| Vee(|r - r'|)|0(r')P , (2) 

where 14e(^) = e^/(er) is the long-range Coulomb 
interaction, with e an effective dielectric constant. 
In the atomic limit the Wannier functions can be 
roughly approximated by a zero-angular-momentum 
wavefunction in the symmetric gauge, (j){r) = 
(27r4)-i/2exp[-rV(44)], where Ib = ^fh^f^ is the 
magnetic length. Simple algebraic manipulations on 
Eq. (2) yield 




implying that, at least asymptotically, U grows propor- 
tionally to \/B. Microscopic details such as the precise 
shape of the confinement potential or the geometry of the 
lattice might affect the result (3) quantitatively but not 
qualitatively: the scaling ex \fB is robust in the asymp- 
totic limit £b 2r, where 2r is the width of the potential 
minima of the artificial lattice [see Fig. la)]. 

The function A{(jj) in a B field is pictorially illustrated 
in red in Fig. 2a). As a comparison, the black dashed line 
labels A{uj) for an unpatterned 2DEG: we distinguish the 
usual Landau level peaks at frequencies ujn = uOc{n^l/2) 
with integer n. In the nanopatterned sample these peaks 
are split into upper and lower Hubbard peaks by strong 
interactions. In this cartoon the measured cyclotron 
mode at cx 5 is an m^er-Landau-level excitation. The 
measured sublinear mode seen in Fig. 1, instead, can 
be neatly explained as an mtra-Landau- level excitation, 
which lies at a frequency cjhb = U /h (x V^, between 
interaction- induced Hubbard peaks [24]. 

Fitting the data labeled by red filled squares in Fig. Id) 
with the simple functional form ficjHB = Q^\/5[T] we 
find a 2 meV, thereby providing a direct measurement 
of the Hubbard- on-site energy scale for our nanopat- 
terned 2DEG. The measured is a factor of two smaller 
than the value extracted from Eq. (3) with the high- 
frequency GaAs dielectric constant e = 13. In Fig. 2b) we 
illustrate a possible two-photon process that contributes 
to the scattering cross section of the HE collective mode. 
The calculated scattering cross section decays exponen- 
tially for sufficiently large values of B (see SM), in agree- 
ment with the data reported in Fig. 3a). In contrast to 




FIG. 3: a) Resonant inelastic light scattering spectra of the Hubbard mode at three values of the magnetic field and T = 1.7 K. 
b) Dependence of the Hubbard mode at T = 1.5 K on the power P (in /xW) of the HeNe laser used to photo-deplete the 
two-dimensional electron gas. From black to blue the electron concentration per site decreases from 8 ± 2 to 3 dz 2 (see SM). 
Data in panel a) are at P = /xW. c) Temperature dependence of the Hubbard mode at B = 5.48 T and P = /iW displaying 
an activated behavior with an activation energy of 0.2 meV. 



the cyclotron mode, indeed, the intensity of the Hubbard 
mode increases up to 5 5.5 T and then collapses ex- 
ponentially at larger fields. 

The Hubbard mode energy exhibits a rather weak de- 
pendence on electron concentration, see Fig. 3b), which 
is decreased by the photo-depletion technique (see SM). 
In the atomic limit the dependence of the Mott-Hubbard 
gap ficjHB on electron concentration is indeed a small ef- 
fect, of the first order in the parameter t/U . In the limit 
of vanishing electron concentration, the strength of the 
transition between the two HBs also vanishes because 
there are no available states in the upper HB (see SM). 
This is in agreement with the large dependence of the 
intensity of the Hubbard mode on electron density re- 
ported in Fig. 3b). Finally, the Hubbard mode displays 
a large sensitivity to temperature changes [see Fig. 3c)] 
and disappears close to T = 5 K. 

We now focus on the low-energy portion of the exci- 
tation spectra, i.e. hw < 1 meV. In ordinary 2DEGs 
this sector is characterized by the spin- wave (SW) mode, 
a spin-flip excitation across the spin gap that, at long 
wavelength, occurs at the bare Zeeman energy ^/xb^t, 
where /iB is the Bohr magneton, g is the Lande gyromag- 
netic factor, and 5t = \J + total magnetic 

field [see Fig. lb)]. The inset to Fig. 4 shows a represen- 
tative result at = 5.5 T. The SW mode is visible at 
energies close to 0.15 meV. The SW energy versus total 
field is reported in Fig. 4 as black filled circles. The spin 
mode is not visible below Bt = 3 T. 

The inset to Fig. 4 displays an additional strong and 
sharp mode above the SW, which has no counterpart in 
an unpatterned 2DEG. The energy dependence of this 
mode is shown in Fig. 4 as red filled triangles. The split- 
ting A of this mode from the SW (filled squares) occurs 
above a threshold B value and depends on the perpen- 



dicular magnetic field only^ a fact that underlines the 
pivotal role of electron-electron interactions. The two 
modes disappear at temperatures approaching T = 1 K. 
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FIG. 4: Energies of the spin- wave mode (black filled circles) 
and the higher energy spin- flip mode (red filled triangles). 
The black dashed line is a linear fit to the data using the stan- 
dard Zeeman formula: ^sw = I^'Imb^t- We find \g\ = 0.42, 
in agreement with the value expected for GaAs. Representa- 
tive examples of the two spin excitation modes at two different 
laser energies (red line: 1522.6 meV; black line: 1522.4 meV) 
are reported in the inset. The black filled squares label the 
splitting A between the two spin modes. 

The observation of a spin doublet is very intriguing 
and suggests the occurrence of a novel correlated state 
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with a gap A. Different types of Coulomb-driven broken- 
symmetry scenarios have been proposed in the context of 
graphene at large magnetic fields [25-29] and linked to 
observations of gap openings in magneto-transport exper- 
iments [13-15]. Remarkably, one of these scenarios [27] 
predicts a splitting of the SW mode similar to what we 
see in our experiment [30] associated to the occurrence 
of lattice-scale order in the honeycomb lattice. 

In the case of graphene, however, the high-field regime 
is not experimentally accessible: one is always in the 
weak-field regime, that is, the magnetic length is much 
larger than the interatomic distance. This is not the 
case in our artificial honeycomb lattice. To support the 
existence of graphene-like effects in our system, we car- 
ried out calculations of the density-of-states based on 
a tight-binding model [31] in the presence of disorder 
comparable to the hopping energy and in the ultra high- 
magnetic-field regime {£b < cl). These results reveal the 
persistence of a structure reminiscent of the zero-energy 
Landau level of graphene (see SM). Similarly to what 
happens in graphene [13-15, 29], electron-electron inter- 
actions can lead to a re-organization of this low-energy 
sector yielding a broken-symmetry ground state with an 
energy gap ^ A. Additionally, the observed opening of 
the gap above a threshold magnetic field indicates a deli- 
cate interplay between hopping, disorder and many-body 
effects. 
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SUPPLEMENTARY MATERIAL 

This section contains technical details and numerical 
results relevant to the main text. 



Nanofabrication of the artificial lattice and resonant 
enhancement of light scattering peaks 

The in-plane potential modulation is achieved by defin- 
ing an array of Nickel disks (with diameter 2r) arranged 
in a honeycomb-lattice geometry (with lattice constant a) 
by e-beam nanolithography and then by etching away the 
material outside the disks by inductive coupled reactive 
ion shallow etching [10, 21]. Owing to the dependence of 
band-bending profiles on GaAs cap-layer thickness, the 



resulting pillars [see Fig. la) in the main text] induce 
a lateral potential modulation with an amplitude of a 
few meV's acting on the electronic system [10]. By tun- 
ing the etching depth d we can reach different regimes 
{i.e. different values of Vq). We focus on a sample with 
d ^ 60 nm for which we estimate Vq ^ 4 meV [11]. 

The resonant inelastic light scattering experiments 
were performed by using a ring-etalon Ti: Sapphire laser 
with a tunable wavelength of circa 800 nm in resonance 
with the magneto-luminescence of the nanostructured 
semiconductor (data not shown) focused on the array 
with a 100 jam diameter area. The scattered light was 
collected into a triple grating spectrometer with CCD de- 
tection: for the geometry of the experiment see Fig. lb) 
in the main text. The intensity of the incident radia- 
tion was kept to values well below 10~^ W/cm^ to avoid 
significant heating of the electrons. 

Crucial for the optical detection of the Hubbard and 
spin modes of the electrons in the honeycomb lattice 
is the resonant enhancement of the corresponding light 
scattering cross sections. This occurs when the incoming 
laser energy matches the energy of an inter-band transi- 
tion of the host crystal. In the experiment this is achieved 
by tuning the laser energy and plotting the light scatter- 
ing signal as a function of energy shift with respect to 
the laser energy. Examples of the resonant effect for the 
Hubbard mode and for the spin- wave and spin-fiip modes 
(reported in Figs. 1, 3, and 4 of the main text) can be 
found in Fig. 5. When displayed as a function of energy 
shift, the collective modes occur at a fixed energy and 
their intensities follow a resonant profile. 

Light scattering spectroscopy of the cyclotron mode 

Representative resonant inelastic light scattering spec- 
tra at T = 1.7 K manifesting the cyclotron mode at two 
values of the perpendicular magnetic field B are reported 
in Fig. 6. The cyclotron mode is parity- forbidden in the 
dipole approximation, and its non-zero albeit weak cross 
section originates from valence-band mixing effects [22]. 
The fact that we see the ordinary 2DEG cyclotron mode 
in our nanopatterned system is due to a magnetic length 
Ib that is smaller than the lattice parameters in the ex- 
plored range of magnetic fields. (The magnetic length 
£^at5 = 2Tis'^18 nm, which is smaller than the 
pillar diameter 2r ^ 60 nm.) There is a second peak 
above Uc (the red solid line in Fig. 6 is a fit with two 
gaussians, which are shown as blue dashed lines) that 
can be associated to a finite- wavevector excitation due to 
the dispersion of the cyclotron mode [22] and the opti- 
cal grating effect of the lattice [32] . The cyclotron mode 
can be detected for 5 ^ 2 T, a value associated with 
electrons occupying the lowest Landau level only as indi- 
cated by the magneto-transport analysis (data not shown 
and Ref. 11). We found that its intensity does not vary 
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FIG. 5: Resonant profile of inelastic light scattering spectra 
showing the Hubbard mode at T = 1.7 K and B = 5.48 T 
[(panel a)] and the spin-wave (SW) and spin-flip (SF) modes 
at T = 55 mK and B = 5.48 T [panel b)]. The incident 
photon energies are indicated. 
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FIG. 6: Representative resonant inelastic light scattering 
spectra showing the cyclotron mode at T = 1.7 K and 
B = 5.48 T and 5.98 T. Incident photon energy is 1547.7 meV 
and 1546.5 meV, respectively. A fit (red solid line) with two 
gaussians (blue dashed lines) is also shown. 



significantly in the explored range of magnetic fields. 



Evaluation of the scattering cross section of the 
Hubbard mode 

In Fig. 2b) of the main text we have shown a possi- 
ble two-photon process that contributes to the scattering 
cross section of the Hubbard collective mode. In this 
scheme the initial state |1) is characterized by a pair of 
neighboring singly-occupied sites with two active elec- 
trons in an antiparallel spin configuration. In the final 
state 1 2) the two initial electrons occupy the same site. 
The final excited state is separated from the ground state 
by the Hubbard charge gap [/, i.e. by the energy cost of 
having two antiparallel spin electrons on the same site. 
Because the lower and upper states split by U have the 
same parity (as they both emerge from the very same 
single-particle level at energy sq) the process depicted in 
Fig. 2B of the main text is allowed in the dipole approx- 
imation, and its intensity is thus expected to be much 
larger than the one of the cyclotron mode at ujc [be- 
cause of the vanishing of the dynamical structure fac- 
tor S{uj) at Uc] in agreement with the experimental ob- 
servation. The scattering cross section can be roughly 
evaluated in the dipole approximation [33] by assuming 
that the valence-band holes are not strongly affected by 
the external periodic modulation. Within the same de- 
gree of accuracy used to derive Eq. (3) in the main text, 
we find that the cross section decays exponentially, i.e. 
d(j/dVt' (X exp [— a^/(4£^)], for sufficiently large values of 
B in agreement with the experimental observation. 



Reducing the electron density by the 
photo-depletion technique 

The electron density of a modulation-doped quantum 
well like the one used in our study can be reduced by 
continuous illumination with photon energy larger than 
the quantum- well (QW) barrier (^ 1.6 eV in our case). 
The mechanism of photo-depletion is based on the fact 
that photoexcited electrons in the AlGaAs barrier con- 
tribute to a charge compensation of the ionized donors 
while photoexcited holes are swept in the QW region, 
thereby reducing the electron density through electron- 
hole radiative recombination [34]. 

In our studies we used a HeNe laser at energy of 1.96 eV 
with powers P up to 50 /iW. The calibration of electron 
density versus HeNe power P is obtained by monitoring 
the evolution of the QW photoluminescence (PL) in an 
unpatterned region of the sample. Representative spectra 
are shown in Fig. 7. The QW optical emission exhibits 
the usual line shape of a modulation-doped QW lumines- 
cence, as determined by the recombination of electrons 
from the bottom of the subband up to the Fermi energy 
(see inset to the left panel of Fig. 7). The distance be- 
tween the main PL peak at energy Ei and the shoulder at 
energy E2 is related to the electron density n through the 
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FIG. 7: (Left panel) Evolution of the photoluminescence spec- 
tra of the unpatterned 2DEG with varying HeNe power P 
(black /iW; red 10 /iW; blue 25 /iW). The in- 
set shows a schematic profile of the conduction and valence 
bands: occupied states are indicated by gray-shaded areas. 
El and E2 refer to the main structures in the photolumines- 
cence. (Right panel) Estimated electron concentration per 
site in the patterned region of the sample versus HeNe laser 
power. 



relation \Ei — E2\ = £^f(1 + ^e/^h)7 where me(h) are the 
effective masses for electrons (holes) and Ep = Trnh'^ /ttIq 
is the Fermi energy. This relation thus allows us to esti- 
mate the electron density n in the unpatterned 2DEG as 
a function of P. This simple analysis cannot be applied 
to the patterned region since in this case the PL line- 
shape is significantly modified as discussed in Ref. 10. 
The 2DEG density n can be significantly reduced from 
its zero-power value ('^ 1 x 10^^ cm~^) as the HeNe power 
P is increased. 

The right panel in Fig. 7 shows the estimated electron 
concentration per site at different HeNe powers in the 
nanopatterned region using the calibration method de- 
scribed above (but starting from the P = value of the 
electron density ~ 3 — 4x 10^^ cm~^ of the nanopatterned 
sample). 



Evaluation of density dependence of the Hubbard 
mode 



In this Section we demonstrate that the Mott-Hubbard 
gap huJuB depends very weakly on the electron concentra- 
tion in the single-band Hubbard model and in the atomic 
limit ([/ > t). 

In the so-called "Hubbard I" approximation [35] , which 
works well at strong coupling, the Hubbard bands E±{k) 
for spin a electrons are given by the following expression: 



E±{k) 



£0 



t{k) + U 
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FIG. 8: Representative inelastic light scattering data showing 
the evolution of the Hubbard mode versus HeNe power at two 
different magnetic fields. The density is decreased (from black 
to blue) by means of the photo-depletion technique. 



where t{k) is the bare band energy and n_cr is the con- 
centration of electrons with spin projection opposite to 
a. In the atomic limit Eq. (4) simplifies to 



\ t{k) 



E^{k)=so^U 
[ E_{k)=eo^t{k){l-n-,) 



(5) 



We thus immediately see that one band is centered at £0 + 
U and the other one at £0, independently of the electron 
concentration. The dependence of the Mott-Hubbard 
gap hwuB = mmk[E^{k) - E_{k)] = U ^ 0{t/U) on 
electron concentration is a small effect, of the first order 
in t/U. In the limit n_cr = the two Hubbard band en- 
ergies are at Sq + t{k) and Eq + [/, but the upper band 
has a vanishing number of states (n_cr) associated with 
it. In other words, the strength of the transition from 
the lower to the upper band vanishes because there are 
no available states in the upper band. 

A similar result is valid also within the multiband gen- 
eralization of the "Hubbard I" approximation [36]. The 
Green's function in the strongly-correlated limit is deter- 
mined by the equation 



G-HE,k) = G-,\E)-t{k) 



(6) 



± ^/t^{k) + U^ +2t{k)U{2r 



1) (4) 



where Gs,t{E) is the energy-dependent Green's function 
of the atomic problem. One can prove (see the Appendix 
in the work by Lebegue et al. [36]) that this expression 
is exact for the multiband Hubbard model up to first 
order mt/U . The poles of this function, which determine 
the positions of the centers of the Hubbard bands, do 
not depend on electron concentration, while the residues, 
which determine the width of the Hubbard bands, do. 
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These theoretical arguments are fuhy consistent with 
the experimental observations shown in Fig. 3b) of the 
main text and in Fig. 8. 



Tight-binding calculations in the 
ultra-high-magnetic-field regime 

We have performed extensive tight-binding calcula- 
tions of the density-of-states (DOS) of non-interacting 
electrons hopping on a honeycomb lattice in the presence 
of a perpendicular magnetic field B and of uniformly- 
distributed disorder of amplitude W. We have simulated 
both on-site disorder, which is mathematically described 
by an Hamiltonian of the type 



n 



dis 



(7) 



with €i uniformly distributed in the interval 
[— W/2,+W/2] and disorder due to random hop- 
ping, which is described by an Hamiltonian of the 
type 



dis 



(id) 



13 3 



(8) 



where 5tij is uniformly distributed in the interval 
[— +W/2]. Details on the numerical technique can 
be found e.g. in Ref. 31. 

In Figs. 9-10 we collect our main findings for a hon- 
eycomb lattice with a = 130 nm and B varying from 
0.001 T up to 5 T. 

In Fig. 9 we clearly see that when is > d the calcu- 
lated DOS is indistinguishable from that of graphene and 
exhibits a zero-energy Landau level and massless Dirac 
fermions at low energies stemming from the topology 
of the lattice [37, 38]. When disorder is switched on, 
the zero-energy Landau level broadens. Fig. 9 illustrates 
the DOS in the ultra-high-magnetic-field regime, i.e. for 
Ib < CL. For B = 2 T and 5 T, for example, we clearly 
see that, in the absence of disorder, a series of Hofstadter 
states appears near zero energy due to commensurability 
effects [39-41] {i.e. the magnetic fiux ^ = 3V3a'^B/2 
through the unit cell of the honeycomb lattice being of 
the same order of the quantum fiux unit hc/e). When dis- 
order is taken into account a broad structure near zero 
energy emerges akin to the zero-energy Landau level in 
disordered graphene (see right panel in Fig. 9). Notice 
that this structure has a width in energy of the order of 
0.2 t for W = 0.5 t, i.e. roughly twice the width of the 
zero-energy Landau level in graphene for the same value 
ofW. 

In analogy to what observed experimentally in 
graphene [13-15, 29], we expect that electron-electron 
interactions can strongly re-organize this low-energy sec- 
tor yielding a correlated ground state with an energy gap 



5 

^ 4 



3 - 
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Q 1 



1 Hill I IIIIH 



1 


1 







-3-2-10 1 2 3 
E/t 



-0.05 0.00 0.05 

E/t 



FIG. 9: Calculated DOS [in units of l/(ta^)] as a function of 
energy E (in units of t) for electrons hopping on a honeycomb 
lattice with lattice constant a = 130 nm and in the presence 
of a small magnetic field B = 0.001 T. We clearly see the 
zero-energy Landau level, unevenly spaced Landau levels of 
massless Dirac fermions at low energies, and equally spaced 
ordinary Landau levels at higher energies. The panel on the 
right is a zoom of the low-energy sector. Black lines label the 
DOS of the clean system, while blue dashed {W = 0.5 t) and 
red dotted (W = 1.0 t) lines label results for the disordered 
system. Data for W have been obtained by using random 
on-site disorder, see Eq. (7). These numerical results have 
been obtained by using periodic boundary conditions on a 
sample containing 3200 x 3200 lattice sites. 



of the order of the spin-doublet energy splitting A in- 
troduced in the main text. Similarly to graphene [29], 
the precise mechanism responsible for the emergence of 
this novel ground state is, to date, unknown. A broken- 
symmetry ground state driven by lattice-scale interac- 
tions beyond the continuum model is compatible with 
the experimental results [27]. 

Before concluding, we would like to emphasize that 
in the limit in which is smaller than a effects be- 
yond those captured by the tight-binding model should 
be taken into account since the shape of the "atomic or- 
bit als" is altered. To a first approximation, the mag- 
netic field simply increases the confinement of electrons 
in each minimum of the periodic potential ("magnetic 
squeezing") thereby renormalizing the hopping ampli- 
tude \t\ 1^(^)1- III the clean limit this fact has no 
implications on the validity of our results in Figs. 9-10. 
Indeed, the data reported in these plots illustrate the 
DOS [in units of l/(ta^)] as a function of E/t: the hop- 
ping amplitude is the only energy scale and its precise 
value thus does not matter. 

In the presence of disorder, however, as magnetic field 
B increases, the tendency to Anderson localization in- 
creases too. The actual phase diagram (and broken- 
symmetry states) of the nanopatterned electron gas re- 
sults from the competition between electron-electron 
interactions, disorder, and external periodic potential. 
This is clearly beyond the scope of the elemental theoret- 
ical analysis reported in the main text and in the SMand 
is left for future investigations. The experimental data 
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FIG. 10: Same as in Fig. 9 but for higher values of the mag- 
netic field, B — 2T (left panels) and B = 5 T (right panels). 
Color coding and labeling are the same as in Fig. 9. Data for 
7^ in the two panels at the bottom have been obtained 
by using disorder due to random hopping, see Eq. (8). 



in Fig. 4 of the original manuscript strongly suggest that 
electron-electron interactions are capable of re-organizing 
the low-energy degrees of freedom of the nanopatterned 
electron liquid at least up to magnetic field values of the 
order of 8-9 Tesla. 
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